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Abstract 

In this paper we consider a new mathematical extension of the Black-Scholes model 
in which the stochastic time and stock share price evolution is described by two 
independent random processes. The parent process is Brownian, and the directing 
process is inverse to the totally skewed, strictly a-stable process. The subordinated 
process represents the Brownian motion indexed by an independent, continuous and 
increasing process. This allows us to introduce the long-term memory effects in the 
classical Black-Scholes model. 
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The option trading has the long history. The mission of the options as finan- 
cial instruments is to protect investors from the stock market randomness. 
Since the early seventies the option market rapidly became very successive 
in development. The theoretical study of options was directed in finding a 
fair and presumably riskless price of these instruments. Without questions, 
the works of Black and Scholes [1] and Merton [2] are a turning-point in the 
study. Their method has been proven to be very useful for investors trading in 
option markets. On the other hand, the approach is fruitful for extending the 
option pricing theory in many ways. Therefore, nowadays the Black-Scholes 
(BS) model is very popular in finance. 

The BS equation is nothing else as a diffusion equation. In fact, their op- 
tion price formula is a solution of the diffusion equation with the initial and 
boundary conditions given by the option contract terms. The fundamental 
principles governing the financial and economical systems are not completely 
uncovered. In recent years the physical community has started applying con- 
cepts and methods of statistical and quantum physics of complex systems to 
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analyze economical problems j3Pf5|6] (and references therein). The improve- 
ment of the BS model itself did not stand still too. As shown in [7], the 
BS equation can be derived using the Stratonovich calculus. The Gaussian 
assumption of the classical BS model based on ideal market conditions sim- 
plifies analytical calculations, but the empirical studies [8] show the effect of 
non-ideal market conditions on the true option price. In particular, the prob- 
ability distribution of returns has heavy tails in contrast to a Gaussian. This 
explains the great interest to various generalizations of classical results. So, in 
[HI the stochastic dynamics of the stock and currency markets is described by 
the fractional Langevin-type stochastic differential equation that differs from 
the standard Langevin equation. The continuous-time random walk (CTRW) 
model is argued to provide a phenomenological description of tick-by-tick dy- 
namics in financial markets [10]. The present paper gives arguments that the 
CTRW model permits ones to generalize the classical BS model. This natural 
extension is based on the general probabilistic formalism of limit theorems. 
The important preference of the approach is its analytical results. We include 
the long-term memory effects in the stochastic process of the BS model. The 
memory effects are characterized only by one parameter. To change it, one 
can control the contribution of memory effects to the model. The classical BS 
model is a particular case of the new model under the complete absence of 
memory. 

The CTRW model is represented by two Markov processes. One of them cor- 
responds to the random waiting-times between successive jumps, another de- 
fines the random space steps. The geometric Brownian motion is a special 
case of the CTRW, where time is deterministic (see below). Let Ti, T2, ... be 
non-negative and independent identically distributed (i.i.d.) random variables 
describing the waiting times between jumps of a walking particle. Assume 
that Ti belongs to the strict domain of attraction of some stable law with 
index < a < 1. This means that there exist 6„ > (ri G N), and the 
sum &„(Ti + ■ — h Tn) converges in distribution to the process having the sta- 
ble distribution with index a. The range < a < 1 is conditioned by the 
support of the time steps Tj on the non-negative semi-axis. In the discrete 
model the internal time r takes on discrete values with an interval 6t such 
that n < [t/6t] < n + 1, where [x] denotes the integer part of x. There 
exists the limit passage from "discrete steps" of the CTRW to "continuous 
steps". The process &[r/5r] E!=/"' Ti under 6t ^ converges in distribution 
to a new process T(r) f ri/"T(l), where f means equal in distribution, and 
T(l) i Ti. The new process is Markovian, strictly a-stable, totally skewed. 
Since T(r) — )• 00 in probability as r — > 00, the sample paths of {T(r)} are 
increasing almost surely (a.s.). The process {T(r)} is self-similar with expo- 
nent i7 = 1/a > 1 HI], i. e. {T(cr)},>o {c^/"T(r)},>o for all c> 0, where 
denotes equality of all finite dimensional distributions. Without loss of 
generality, we may assume jumps in the one- dimensional space. Denote by Ri 
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the space steps. Let Ri, R2, ... be i.i.d. random variables independent of {Tj} 
and have the Gaussian distribution. Using the hmit passage from "discrete" 
to "continuous" jumps, we obtain the stochastic process {-R(r)}^>o with the 
self-similar relation {R{ct)}t>o ^f' {c^^'^R{t)}t->q for any c > 0. It should 
be pointed out that both processes {-R(r)}^>o and {T{T)}r>o depend on the 
continuous internal parameter r that differs from the real time t. 

To build the continuous position vector of the walking particle, we need the 
process which represents the continuous limit of the discrete counting process 
{Nt}t>o- For t > the number of jumps up to time t is Nt = maxjn G N | 
J2i=i Ti < t}, and the vector r^r^ = X^i^i Ri defines the position of the particle 
at time t. It turns out that the scaling limit of {Nt}t>o is the hitting process of 
{T(x)}x>o- The hitting time process is well defined S{t) = mf{x \ T(x) > t} 
and depends on the true time t. The two processes {T{x)} and {5'(t)} are the 
inverse of each other, S(T{t)) = r a.s. Since {T{x)}x>o is strictly increasing, 
the process {S{t)}t^Q is non-decreasing. From the self-similarity of {T{x)} it 
follows the same property for {S'(t)}, i. e. {S{ct)}t>o {c°'S{t)}t>o for any 
c > 0. While {T(x)}s>o is a Levy process, the inverse process {5'(t)}f>o is 
no longer a Levy process, neither a Markov process, but it is a continuous 
submartingal, as shown in [T2]. The random value S{t) has a Mittag-Leffler 
distribution with (e-^'^W) = E^=o(-^^")"/r(l + = ^a(-fr), where {X) 
denotes the expectation of a real valued random variable X, and T[z) is the 
Gamma function. The sample paths of {Nt}t>o and {5'(t)}(>o are increasing. 
Then the position of the particle at the given real time t is defined by 
the subordinated process R{S{t)). Recall briefly that a subordinated process 
Y{U{t)) is obtained by randomizing the time clock of a random process Y{t) 
using a new clock U{t), where U{t) is a random process with nonnegative 
independent increments. The resulting process Y{U{t)) is said to be subor- 
dinated to called the parent process, and is directed by U{t), called 
the directing process. The directing process is often referred to as the ran- 
domized time or operational time [13]. In general, the subordinated process 
Y{U{t)) can become non-Markovian, though its parent process is Markovian. 
The process R{S(t)) is self-similar with index a/2 such that {R{S{ct))}t>o 
f-d. jc"/^ i?(S'(t))}t>o is for all c > 0. In fact, the position vector = I3s(t) 
represents the randomization of the internal time r of a Brownian motion Br 
by an independent, positive and non- decreasing process S{t). 

The probability density of the position vector with t > satisfies 

00 

/'(t,x)= Jp''{T,x)p'{t,T)dT, (1) 


where p^{t,x) represents the probability to find the parent process R{t) at 
X on the operational time r, and p^{t, r) is the probability to be at the op- 
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erational time r on the real time t. The Laplace transform of the probability 
density of the random variable S{t) with respect to x gives 

oo 

p^it.v) = J e-''^p^{t,x)dx = (e-"^(*)) = Ea{-vr). 





We need also the Laplace transform of p^{t,x) with respect to t. The Mittag- 
Leffler function Eal—vt"') has the following Laplace transform u°'~^/{u°' + v) 
with respect to t. To invert the latter analytically, we obtain 

oo 

p^{u,x) = J e'''^p^{t,x)dt = M°"^exp{-M°a;}. 





In Laplace space the probability density p^^{t,x) has the most simple form 
u°'~^p^{u°',x), as p^{u°',x) = J^P^{t,x) exp{—u°'T} dr. For our purpose, it 
is useful to find the explicit form of the probability density p^^{t, x). According 
to the inverse formula applied to p^{u, x), we have 



p^{t,x) 



1 



(2) 



Br 



where Br denotes the Bromwich path. Make the variable transform ut ^ u 
and denote w = x/t". Then we deform the Bromwich path into the Hankel 
path Ha for which a contour begins at m = — oo — m (a > 0), encircles the 
branch cut that lies along the negative real axis and comes to the end at 
u = — oo + ib {b > 0). Expanding function exp{— w-u"} in a Taylor series 
about w and using the Hankel representation of the reciprocal of the Gamma 
function 

= — [ e'^u'^'du, 
T{z) 2Txi J 

Ha 

we get the following series 

Further, we briefly consider the character of Fa{z). 

The function Fa{z) is an entire function in z. It has the H-function repre- 
sentation (z I ''^o"')"'') [S]- The important property of Fa{z) is that it is 
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non-negative for 2; > 0. It is easily verified that Fa{z) dz = 1. Thus, the 
function can be a probabihty density. The case a = 1 corresponds to the Dirac 
(5-function, Fi{z) = 6{z — 1). In particular cases a = 1/2 and a = 1/3 we have 
Fi/2iz) = exp{— z'^ / A} / and -^1/3(2;) = 3^/^Ai(2;/3^/^) respectively, where 
Ai denotes the Airy function [15]. The function Fa{z) has also other interest- 
ing properties. For < a < 1/2 the function is monotonic decreasing, whereas 
for l/2<Q;<lit has a maximum value at a certain point Zj^ax depending on 
a. It should be observed here that the basic Cauchy and Signaling problems 
of the time fractional diffusion- wave equation can be expressed in terms of the 
function [MTTj . 

Turning back to Eq. ([1]), the probability density p^^{t,x) is written as 

oo 00 

J V nJJt J y z 



where D is the constant. This function is non- negative and satisfies the nor- 
malization condition 



x) dx = J Fa{z) dz = 1. 



Since the parent process {R{t)} and the directing process {S{t)} have finite 
moments of any order, the subordinated process {R{S{t))} has finite moments 
of any order too. The first and second moments of can be obtained by the 
direct calculations: 

(r*) = , 

00 „ 
1 r nt" 

(r^) = -Dr z FJz) dz = . 

^ *^ 2 } ^ ' 2r 1 + a 

^ ' 

The process behaves as subdiffusion (0 < a < 1). Note that the boundary 
case a = 1 may be also included in the consideration because of T(r) = r 
a.s. Then the hitting time process is deterministic, S{t) = t. The probability 
density p^{T,t) degenerates in the Dirac (5-function so that x) becomes 
equal to p^{t,x). The constant D is interpreted as a generalized diffusion 
coefficient with dimension [D] = length^/ time". 

The ordinary Brownian motion satisfies the stochastic differential equation 
(SDE) 

dR{T) = f{R{T)) dr + g{R{T)) dBr , 
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where / and g are some functions. The process subordinated to the Brown- 
ian motion {R{S(t))}t>o = ^s{t) is a continuous martingal and the directing 
process {5'(t)} is a continuous submartingal with respect to an appropriate 
filtration [IB]. Therefore, the subordinated process obeys the following SDE 

dvt = f{vt)dS{t) + g{rt) dBsit). 



In the classical BS model the evolution of the option price is governed by the 
Brownian motion. The well-known BS formula is of the form 

C(r, x) = x^d+) - Ke-^^^{d_) , /3 = 2r/a'^ , 



where x is the share price, K the striking price, r the interest rate, a the 
volatility, and the probability integral 

1 r 

Hz) = ^ / exp{-yy2}dy 

— oo 



is calculated for 

d± = {2r)-^^^\\n{x/K) +t{P±1) 



If the price evolution is consequent of the subordinated process Bs(t), the BS 
formula transforms into 

oo 

S{t, x) = t"" J F^z/t") C{z, x) dz . (4) 



At a = 1 we obtain the classical BS formula. All financial derivatives (op- 
tions of any kind, futures, forwards, etc.) have the same boundary conditions, 
but different either initial or final condition [19]. The detailed comparison of 
the various cases for this new model (HI) will carry out elsewhere. The frac- 
tional extension of the BS model has been considered also in [20], but on the 
macroscopic basis without any microscopic dynamics presented above. 

Finally, we note that the index a characterizes memory effects in the sub- 
ordinated process r^. Let L{x) be a time-independent Fokker-Plank opera- 
tor, whose exact form is not important for the following. If the ordinary 
Fokker-Plank equation (FPE) dp^{T,x)/dT = [L{x) p^]{t,x) describes the 
evolution of a Brownian particle, the probability density p^^{t,x) satisfies the 
fractional FPE. This can be shown by simple computations. Using the rela- 
tion p^^{u,x) = u"'~^p^{u"',x) in Laplace space and acting the operator L{x) 
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on the Laplace image [L(x) p'"*]('U, x) takes the form u°'p^*{u,x) — 

f{x)u°'~^, where f{x) is the initial condition. The inverse Laplace transform 
of the latter expression gives the above-mentioned fractional FPE 

1 I 

= f{x) + —-j dr{t - r)"-i[L(x)p^*](r,x) . 
^ ' 

The kernel of this integral equation is a power function. It just causes the 
long-term memory effects in the process of interest. As shown in [l7j, due to 
such kind of memory effects, the complex nature of the microscopic behavior 
of stochastic systems can be transmitted to the macroscopic level of their 
dynamics. 
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